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1.  Introduction.  A  chemical  vapor  deposition  (CVD)  process  has  been  used  to  coat  short 
length  single  crystal  sapphire  fibers  with  P-AI2O3.  In  order  to  make  the  process  viable  for 
producing  the  longer  coated  fibers  needed  for  practical  composite  materials,  it  will  be 
necessary  to  use  a  continuous  CVD  process.  This  will  require  the  determination  of  critical 
parameters  and  a  control  strategy  for  optimal  coating. 

We  have  developed  a  mathematical  model  of  the  coating  process  and  used  our  adap¬ 
tive  finite  element  software  to  perform  a  parameter  study  of  film  coatings  and  their  growth 
rates  as  functions  of  operating  temperature,  precursor  mass  fraction  and  distribution,  sys¬ 
tem  pressure,  flow  rate,  fiber  speed,  reaction  probability,  and  reactor  geometry.  This  ini¬ 
tial  study  has  helped  to  determine  critical  parameters  of  the  process  and  will  help  to 
develop  production  software. 

The  system  consists  of  a  carrier  gas  entering  at  the  bottom  of  a  cylindrical  hot-wall 
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reactor  containing  a  single  central  fiber  and  exiting  at  the  top.  The  fiber  may  be  drawn  in 
either  the  same  or  the  opposite  direction  as  the  gas  flow.  Both  the  gas  mixture  and  the 
fiber  enter  the  reactor  at  ambient  temperature.  When  the  flow  stablizes,  small  quantities  of 
the  precursor  are  added  to  the  carrier  gas.  The  gas  is  heated  by  a  combination  of  conduc¬ 
tion  and  convection  while  the  fiber  is  heated  by  convection  and  radiation.  A  thin  layer  of 
P-AI2O3  coats  all  heated  solid  surfaces.  Typical  values  of  operating  parameters  and  reac¬ 
tor  geometry  are  given  in  Table  I. 

Table  I.  Typical  Operating  Parameters  and  Reactor  Geometry. 


Reactor  diameter 

2.5  cm 

Reactor  length 

1  m 

Fiber  diameter 

40  pw 

Oxygen  flow  rate 

250  5cc/niin 

Argon  flow  rate 

100  see  1  min 

Outlet  pressure 

10-20  torr 

Reactor  Temperature 

900°C 

2.  Mathematical  Model.  We  model  the  system  as  a  three-dimensional  axisymmetric 
steady  flow  of  an  ideal  gas.  As  an  initial  step,  we  neglect  reactions  in  the  fluid,  the  effect 
of  the  deposited  films  on  the  fluid  flow,  thermal  diffusion,  and  the  heat  released  by  surface 
reactions.  We  further  assume  that  every  precursor  molecule  that  collides  with  the  wall  or 
the  fiber  surface  has  a  probability  a*  of  reacting  and  depositing  A/2O3  on  the  surface  and 
that  inert  products  are  carried  away  by  the  flow.  At  any  instant,  the  mixture  contains  only 
oxygen,  argon,  and  precursor  which  are  indexed  by  i  =  1,  2,  3,  respectively.  The  deposited 
A/2O3  is  indexed  by  i  =4.  Let  and  x,  ,  respectively,  denote  the  mass  and  mole  frac¬ 
tions  of  species  i  =  1,  2,  3,  present  in  the  gas  mixture  and  Af,  denote  the  molecular 
weight  of  species  i  =  1,  2,  3,  4. 

2.1.  Fluid  Flow.  The  mathematical  model  [1,2]  of  the  gaseous  mixture  in  cylindrical  coor¬ 
dinates  (r,  z),  consists  of 


*  A  list  of  symbols  appears  in  Table  11  near  the  end  of  this  report 
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(i)  a  continuity  equation 


(ii)  momentum  equations 


1  ^(^pVr)  ^(pVz) 
r  dr  ^  dz 


P^r-^  +  ^z 


dv^  _  dp  1  ^(^V) 


dz  ' 


(2.2a) 


dp  1 

dz'^  r  dr  dz 


(iii)  stress  tensor  components 


(2.2b) 


2 


dr  r  dz  '  3  ^  dr  ^  r  9z  ’ 


(2.2c,d) 


Vz  =  P 


(iv)  an  energy  equation 


3v^  dv^ 


2 

=  1“  'it 


3T  dT  1  a  ,  ar  a  ,  ar 
^'><’'’'-3;^^'’^^  =  71?*^  *Tz'‘lI  ' 


(v)  a  reacting  species  transport  equation 


(2.2e,f) 


dY^  dY^ 

and  (vi)  an  equation  of  state 


13  3  ^^3 

7¥  ^  ’ 


p=^. 

RJ 


(2.5b) 


The  quantities  and  are  the  radial  and  axial  components  of  the  mixture  velocity  v,  p 
is  the  pressure,  T  is  the  temperature,  p  is  mixture  density,  p  is  the  mixture  viscosity,  k  is 
the  mixture  thermal  conductivity,  g  is  gravitational  acceleration,  Cp  is  the  mixture  specific 
heat,  D  is  the  precursor  diffusivity  in  the  mixture,  Rg  is  the  universal  gas  constant,  and  M 
is  the  mixture  mean  molecular  weight  given  by 


M  = 

i=l 


(2.5b) 


assuming  X2«  Xi,  i  =  1,  2. 


2.2.  Surface  Reaction  and  Film  Growth.  As  a  first  approximation  in  modeling  the  com¬ 
plex  surface  reaction,  we  have  assumed  that  every  precursor  molecule  that  hits  the  surface 
reacts  according  to 


2{Al-{0-C2H2)'i)gas  +  3(02)ga5 - ^  ^(^20)gas^ 

(2.6a) 

to  deposit  the  desired  material.  By-products  of  this  reaction  are  assumed  to  be  inert  and 
carried  away  by  the  flow.  The  impingement  rate  I  of  the  precursor  in  the  mixture  is  the 
number  of  molecules  of  this  species  that  hit  the  surface  per  unit  area  per  unit  time,  i.e.. 


I  = 


PX3 


(2.6b) 


yjlnm^KT 

where  k  is  Boltzmann’s  constant  and  m3  is  the  mass  per  molecule  of  the  precursor. 


Assuming  the  flux  of  precursor  reacting  on  the  solid  surfaces  to  be  proportional  to  the 


impingement  rate  I  and  the  reaction  probability  ol(T),  and  using  K  =  Rg/N^y  and  the  mix¬ 
ture  relation  [6] 


Y3/M2 

with  j  =  1,  2,  we  obtain 


where 


Y3P 

/,  =pa(r)-^ 


P  = 


{M/M2) 


aiT)  =  Ae 


^2TtM2Rg 

and  the  activation  energy  AH  and  A  are  to  be  selected  (see  Section 


(2.6c) 


(2.7a) 


(2.7b,c) 


We  denote  the  thickness  of  the  film  on  the  fiber  surface  by  h{z)  and  use  (2.7)  to 


obtain 
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where  P4  is  the  density  of  A/2O3  and  L  is  the  reactor  length. 

2.3.  Boundary  Conditions.  The  system  (2.1-8)  is  closed  by  adding  boundary  conditions  at 
the  inlet,  outlet,  wall,  and  fiber  surface.  We  consider  cases  involving  (i)  flow  in  the  whole 
reactor  and  (ii)  flow  in  the  central  part  of  the  furnace. 

2.3.1.  Whole  Reactor  Boundary  Conditions.  To  determine  the  flow  in  the  entire  tube,  we 
solve  (2.1-8)  subject  to  the  following  conditions: 


(i)  at  the  inlet 


(ii)  on  the  wall 


Vr  =0,  V,  =  V;„,  T  =  Ti„,  73  =  p  =Pi„, 
at  z  =  0,  Rj^  <  r  <  R, 


(2.9a,b,c,d,e) 


9^3 

=  0,  V  =  0,  T  =  T^iz),  =  M^I^,  at  r  =R,0<z  <L, 

dr 

(2.9f,g,h,i) 


(iii)  at  the  outlet 


dz 


=  0, 


973 

9z 


=  0,  at  z  =  L,  Ry  <  r  <  R,  (2.9j,k,l,m) 


and  (iv)  on  the  fiber 


=  0,  v^=Vf,  pD  =M7,Ir,  h{Qi)  =  0,  (2.9n,o,p,q) 

p/  '=pf  -  /c  =  O'  ^<«/'  o>  =  =  O' 

at  r=Rf,0<z^L,  (2.9r,s,t) 


where  radiant  heat  flux  from  the  wall  is  neglected,  =  kdTIdr  is  the  convective  heat 
flux  from  the  gas  to  the  fiber,  and  Rjr  and  Ry^,  are  the  fiber  and  reactor  radii,  respectively. 
Heat  conduction  in  the  fiber  is  modeled  by  the  one-dimensional  heat  conduction  equation 
(2.9r,s,t)  which  is  coupled  to  the  gas  energy  equation  through  /^.  The  quantities 
p/>  Cpjr,  kf,  and  Vy  denote  the  fiber  density,  specific  heat,  thermal  conductivity,  and  velo- 
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city,  respectively,  while  v,„  r,„ ,  ,  w,„  are  the  axial  velocity,  the  temperature,  the  pres¬ 

sure,  and  the  precursor  mass  fraction  at  the  inlet  of  the  reactor,  respectively. 

2.3.2.  Central  Core  Boundary  Conditions.  Assuming  the  flow  to  reach  a  constant  tem¬ 
perature  T  away  from  the  inlet  and  outlet,  the  density  becomes  constant  in  the  central  part 
of  the  furnace.  Therefore,  we  model  the  central  part  of  the  CVD  reactor  by  (2.1-2),  (2.4-8) 
subject  to  the  boundary  conditions: 


3^3 

h(Zi„)  =  0,  Vr  =  0,  v^=Vy.,  pD 

at  r  =  Rf ,  0  <  z  ^  L.  (2.10k,l,m,n) 

2.4  Physical  Properties.  Evaluation  of  physical  properties  of  the  gas  mixture  is  per¬ 
formed  using  experimental  tables  and  mixture  theory  [6]  assuming  an  argon  and  oxygen 
mixture,  i.e.,  X3  =  0.  The  mixture  viscosity,  thermal  diffiisivity,  and  specific  heat  are  given 
by 


il  =  X  >  ^  =  i  2'^'  ’  (2.11a,b,c) 

;=i  y=i 


with 
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<t>n-  = 


r  11/2f  '|l/4‘]2 

1+  Mj/Mj 

_ J _ _ J 

V8[l  + 


(2.1  Id) 


where  jO.,- ,  ,  and  c^,  denote  the  viscosity,  thermal  conductivity,  and  specific  heat,  respec¬ 

tively,  of  species  /  =1,  2. 

Using  a  least-squares  fit  to  the  data  generated  from  (2.11)  and  the  experimental  tables 
[6],  we  obtain 


^(r)  =  10"^  (0.351796  -h  0.68085210"^  T  -  0.301 185510"'’  -h  0.818010710"’' 


Cp{T)  =  103(0.727147  +  0.19714910-3 T  -  0.654446410-^73  _  o.26141910-*Or3), 

k{T)  =  10-3 (0.308270  -h  0.75678010"®  T  -  0.139484010-9  T^).  (2.12a,b,c) 

The  diffusion  coefficient  of  precursor  in  a  mixture  of  oxygen  and  argon  is  computed 
using  the  binary  diffusion  theory  [6] 


D=m,inr^-^ 


3,2  Na,  ^|T\l/Mi+l/M) 


(2.1 2d) 


We  estimate  the  collision  diameter  5  =  12.5  10  to  obtain 


D  =  cj - ,  with  Cj  =  2.5610  ^ 

P 


(2.12e) 


3.  Numerical  Solution  Procedure..  We  obtain  a  dimensionless  version  of  (2.1-10)  by 
using 


r  =  Lf,  z  =  Lz,T  =  Tf,  =  Uv^,  =  Uv^,  T3  =  Wy3,  p  =  pUy, 


h  =  Hh,  k  =  kk,  D  =  DD ,  Cp  =  CpCp,  and  |X  =  |ijl. 


(3.1a) 


where 


r  =  r,.„,  L  =  L,  W  =  1,  p  =  H  =  \m,k=  k(Ti„),  Cp  =  Cp(Ti„), 


P  =  ),  and  D  =  D  (7)„ ,  )  (3.1b) 

In  the  remainder  of  this  report,  we  consider  the  dimensionless  system  and  omit  the  ~. 
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3.1  Stabilization  Method.  Since  standard  Galerkin  finite-element  methods  based  on  equal 
polynomial  degrees  for  both  velocity  and  pressure  fail  to  satisfy  the  Babuska-Brezzi  stabil¬ 
ity  condition  [4],  we  implemented  the  Douglas-Wang  stabilized  streamline  upwind  method 
[5]  that  consists  of  perturbing  the  standard  finite  element  formulation  by  the  elemental 
least-squares  term 


A^a 


=  E  ,  Phi  Ck i\h )  A[V;, ]iqh ,  w,, ) 

K  =  1^ 

where  the  L2  inner  product  on  element  K  is 


(3.2a) 


((t),  =  l^VfdK, 

K 

\f^  and  Wfj  denote  trial  and  test  functions,  respectively,  for  the  velocity;  pf^  and  denote 
the  trial  and  test  functions,  respectively,  for  the  pressure;  R[v,  p  ]  is  the  residual 


R[v,p]  = 


dVy 


dz 

dv. 


dv, 


dz 


dp _ 1_ 

dr  Re 

dp _ 1_ 

dz  Re 


1  djrx^)  ^  ^ 
r  dr 


dx 


r  dz 


rz 


1  ^'^zz 


dr 


dz 


^fr 


,  (3.2b) 


where  Re  =  pC/L/|i  is  Reynolds  number,  and  Fr  =  U^/{gL)  is  Proud  number.  Finally,  the 
operator  A[v;,  ]  is  given  by 


A[v](^,  w)  = 


dwj. 

^  -f-  - 


dr 


dz 


9w, 


dz 


dq  1 

dq  1 


1  d{r^^)  % 


d^ 


rz 


r  dr 
1 


d\ 


dz 


zz 


dr 


dz 


(3.2c) 


where  and  ^  are  obtained  from  equations  (2.2c-f)  with  v  replaced  by 

and  Y  =  —  1.  There  are  two  other  popular  choices  of  y:  the  choice  y  =  0  yields  the  Stream¬ 
line  Upwind  Petrov-Galerkin  (SUPG)  method  [7,8]  and  7  =  1  yields  the  Galerkin-Least- 
Squares  method  [9]. 


Letting  hfr  be  the  diameter  of  element  K  and  following  Franca  and  Frey  [7]  yields 


(3.2d,e) 
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where 


Hk 

Cj^(v)  =  0.01  with  XfT  = 


Ivlj  =  (V,2  +  REi^  = 


RE  \y\2hic 
12 


,  and  ^(REj()  = 


jREf^,  if  0  <  RE^  <  1 

|l,if/?£j^  >  1. 


(3.1f,g,h) 


3.2  Adaptive  Software.  The  resulting  system  modeling  the  reaction,  convection,  and 
diffusion  occurring  during  the  CVD  processing  is  solved  using  adaptive  software  that  can 
reduce  computational  difficulties  while  automating  many  of  the  decisions  associated  with  a 
numerical  solution.  Our  adaptive  system  has  capabilities  for  automatic  quadtree-structured 
mesh  generation,  mesh  refinement/coarsening  (h-refinement),  method  order  variation  (p- 
refinement),  and  mesh  motion  (r-refinement);  however,  only  h-refinement  was  used  in  this 
initial  study.  The  mesh  is  organized  in  a  tree  structure  where  the  root  node  is  the  domain 
having  as  offspring  all  elements  of  a  base  mesh.  An  error  indicator  is  associated  with  each 
quadrilateral  element  and  h-refinement  is  used  to  refine  those  elements  having  large  errors. 
Elements  may  be  divided  into  four  sub-elements  and  these  are  regarded  as  offspring  of 
their  coarser  parents  in  the  tree  structure.  Leaf  nodes  of  the  tree  represent  the  actual  ele¬ 
ments  of  the  mesh.  Solutions  of  two-dimensional  stationary  and  transient  problems  begin 
with  a  trial  solution  computed  on  a  coarse  mesh  using  a  first-order  method.  The  solution  is 
halted  periodically  and  an  estimate  of  the  discretization  error  is  appraised.  Meshes  are 
subdivided  in  high-error  regions  in  an  attempt  to  maintain  the  error  estimate  within 
prescribed  limits.  Estimates  of  spatial  discretization  errors  are  obtained  from  jumps  in 
gradients  of  the  finite  element  solution  at  inter-element  boundaries  [3]. 


4.  Computational  Results 

The  adaptive  software  system  described  in  Section  3  has  been  used  to  solve  a  dimension¬ 
less  version  of  (2.1-10)  using  a  piecewise  bilinear  finite  element  solution  on  quadrilateral 
meshes.  The  initial  base  mesh  was  graded  near  the  fiber  and  the  wall  to  improve  the 
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convergence  rate  in  boundary  layers.  A  25x25  base  mesh  was  used  with  up  to  four  levels 
of  refinement.  The  resulting  discrete  systems  were  solved  using  a  quasi-Newton  iteration 
procedure  for  steady  systems  or  the  backward-difference  software  DASSL  for  transient 
systems  [10].  In  all  cases,  we  selected  Tj  =  0.667,  Y2  =  0.333,  and  40|i  diameter  fibers 
drawn  at  10~^ m /sec. 

4.1.  Nonisothermal  Case.  First,  we  present  results  for  the  whole  furnace  using  (2.1-9). 

We  chose  L  =  30cm,  R  =  1cm,  p  =  10 torr,  dev  =  0.3,  w  =  0.05,  =  300Ar, 

=  nooK, 

T^(z)  =  Ti^+if^-Ti„)s(z)  (Ala) 

with 

5(^)  =  1{1  +  tanh[50(z  -  0.2)]}  {1  -  tanh[50(z  -  0.05)]},  (4.1b) 

and  a  reaction  probability 

AH  AH 

a(T)  =  0.02  e  e  ^  ,  with  AH  =  10\  (4.2) 

The  resulting  pressure,  axial  velocity,  temperature,  and  precursor  mass  fraction  are 
shown  in  Figure  1.  The  gas  and  the  fiber  reach  the  wall  temperature  in  the  central  part  of 
the  furnace.  As  the  gas  is  heated  it  accelerates  in  the  central  part  of  the  furnace.  We  also 
present  the  fiber  coating  thickness  and  deposition  rates  on  the  wall  as  function  of  z  in  Fig¬ 
ure  2.  In  the  preheat  zone  near  the  inlet,  little  if  any  coating  results;  however,  as  soon  as 
the  gas  reaches  the  heated  zone,  surface  reactions  are  triggered  and  the  deposition  rate 
increases  to  reach  a  maximum  and  subsequently  decays  as  the  precursor  is  depleted.  As  a 
result  of  these  computations,  it  appears  that  the  exponential  dependence  of  a  on  tempera¬ 
ture  leads  to  negligible  surface  reactions  for  T  <  T.  We  thus  decided  to  simplify  this  ini¬ 
tial  study  and  confine  our  attention  to  the  central  part  of  the  reactor. 

4.2.  Isothermal  Case.  We  consider  the  central  part  of  the  furnace  and  perform  a 
parametric  study  by  solving  the  system  (2.1-2),  (2.4-8),  (2.10)  assuming  the  wall,  fiber. 
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and  gas  to  be  at  uniform  temperature.  We  consider  cases  where  the  precursor  is  injected 
near  the  fiber  surface.  Let  Fy  =  wy//n,„,  and  F^  =  fn^lnij^  denote  the  per¬ 

centage  of  precursor  that  deposits  on  the  fiber,  deposits  on  the  wall,  and  leaves  the  furnace 
unreacted,  respectively.  Here,  denotes  the  precursor  mass  flux  entering  the  furnace, 
my  the  precursor  mass  flux  depositing  on  the  fiber,  the  precursor  mass  flux  depositing 
on  reactor  walls,  and  m^  the  unreacted  precursor  mass  flux  exiting  the  furnace.  We  use  the 
percentage  Fy  of  precursor  that  deposits  on  the  fiber  to  measure  the  efficiency  of  the  sys¬ 
tem.  The  production  rate  is  measured  by  the  fiber  coating  thickness  at  the  outlet.  At  the 
inlet,  the  precursor  is  injected  within  a  distance  dev*{R-Rj:)  of  the  fiber  surface  and  Wj^ 
is  given  by 

=  ^{1  -  tanh[e(r  -  (Fy  -i-  {R^  -  Rf)dev))]},  with  e  =  .  (4.3) 

l  -  Rf 

Thus,  when  dev  =  1,  the  precursor  fraction  is  uniform  at  the  inlet,  and  when  dev  =  0.5  the 
precursor  is  concentrated  within  a  distance  (R-Rf)/2  of  the  fiber. 

In  order  to  examine  the  effect  of  flow-rate  variations  on  efficiency  and  production 
rate  of  the  coating  process,  we  solved  (2.1-2),  (2.4-8),  (2.10)  for  p  =20torr,  a  =  0.02, 
R  =  1cm,  L  =  20cm,  w  =  0.05,  T  =  1200^,  dev  =  0.05,  and  flow  rates  of  61,  124, 
188,  251,  315,  442scclmin.  The  results  are  shown  in  Figure  3.  Higher  flow  rates  yield 
lower  efficiency  and  higher  production  rates. 

In  order  to  examine  the  effect  of  precursor  distribution,  we  solve  the  system  (2.1-2), 
(2.4-8),  (2.10)  with  p  =  20torr,  a  =  0.05,  F  =  1cm,  L  =  20cm,  w  =  0.05,  T  =  1200F, 
and  dev  =  0.05,  0.1,  0.2,  0.3,  0.4 , 0.5 , 0.6,  0.9,  with  a  flow  rate  of  442scclmin .  The 
results,  shown  in  Figure  4,  indicate  that  the  closer  the  precursor  to  the  fiber  the  more 
efficient  the  system.  This  is  expected  because  the  precursor  has  to  diffuse  further  to  reach 
the  wall.  However,  smaller  values  of  dev  with  a  constant  w  and  flow  rate,  lead  to  smaller 
mass  fluxes  of  precursor  at  the  inlet;  thus,  slower  production  rates  arise  (see  Figure  4, 
left). 
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In  order  to  examine  the  effect  of  reactor  length,  we  solved  the  system  (2.1-2),  (2.4-8), 
(2.10)  for  p  =  lOtorr,  a  =  0.02,  R  =  1cm,  w  =  0.05,  dev  =  0.1,  T  =  HOOK,  a  flow  rate 
of  MAscdmin,  and  L  =  10,  20,  30,  40,  50cm.  The  results  presented  in  Figure  5  show 
that  longer  reactors  improve  both  efficiency  and  production  rate.  However,  we  expect  that 
for  a  given  flow  rate,  temperature,  pressure,  and  reaction  probability,  we  can  determine 
the  optimal  reactor  length. 

In  order  to  examine  the  effect  of  reaction  probability  a,  we  set  p  =  lOtorr , 
R  =  1cm,  L  =  20cm,  dev  =0.2,  w  =0.05,  T  =  1200^r,  a  flow  rate  of  llAscdmin, 
with  a  =  0.001,  0.01,  0.05,  0.1  and  solve  the  system  (2.1-2),  (2.4-8),  (2.10).  We  plot 
Pf ,  and  Pg  vs.  reaction  probability  a  and  the  coating  thickness  /z  as  a  function  of  z  in 
Figure  6.  These  results  indicate  that  for  large  reaction  probabilities  the  precursor  reacts 
quickly  on  the  fiber  surface  as  it  enters  the  reactor  and,  therefore,  only  a  small  percentage 
of  it  reaches  the  wall  or  is  convected  out  of  the  reactor.  Thus,  a  higher  reaction  probability 
increases  efficiency  and  production  rate. 

In  order  to  examine  the  effects  of  precursor  fraction,  we  solve  the  system  (2.1-2), 
(2.4-8),  (2.10)  for  p  =  lOtorr,  a  =  0.02,  R  =  1cm,  L  =  20cm,  dev  =  0.2,  T  =  1200.K', 
and  a  flow  rate  of  633scdmin,  with  w  =  0.05,  0.1,  0.2,  0.5.  The  computational  results 
predict  that  efficiency  is  independent  of  precursor  mass  fraction  for  a  fixed  distribution 
(see  Figure  7,  left).  However,  production  rate  depends  linearly  on  precursor  mass  fraction. 

In  order  to  examine  the  effect  of  reactor  radius,  we  assumed  the  precursor  to  be 
injected  within  a  fixed  distance  from  the  fiber  for  all  values  of  and  chose 
dev  =  0.02(0.5  -  Rf)/(R^^,  -  Rj).  The  other  parameters  are  p=20torr, 
a  =  0.02,  w  =  0.05,  L  =  20cm,  T  =  HOOK,  a  flow  rate  of  350scdmin  and  R^  =  0.5, 
1,  1.5.  The  results,  shown  in  Figure  8,  suggest  that  efficiency  and  production  rate  of  this 
CVD  process  are  improved  by  making  the  reactor  wall  farther  away  from  the  fiber.  The 
flow  velocity  decreases  with  increasing  R^^,,  when  holding  dev  and  w  fixed.  This  leads  to 
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a  lower  precursor  flux  entering  the  reactor.  Consequently,  the  production  rate  of  the  pro¬ 
cess  deteriorates  (see  Figure  7,  right).  Varying  the  reactor  radius  produces  the  same  effect 
as  varying  dev . 

In  order  to  examine  the  effect  of  temperature,  we  assumed  a  temperature-independent 
reaction  probability  a  =  0.02  and  selected  p  =  \0torr,  R  =  1cm,  L  =  20cm,  dev  =  0.3,  a 
flow  rate  of  350scc Imin ,  and  T  =  800,  900,  1000,  1100,  1200,  1300Jfsr.  Results,  shown  in 
Figure  9  (left),  suggest  that  efficiency  is  temperature  independent.  Higher  temperatures 
help  precursor  molecules  diffuse  faster  towards  the  walls  and  deposit  before  depositing  on 
the  fiber  or  exiting  the  furnace  (see  Figure  9,  left).  The  results,  shown  in  Figure  9  (right), 
indicate  that  production  rates  decrease  with  increasing  temperature. 

We  examine  the  effect  of  pressure  on  the  CVD  process  by  solving  the  partial 
differential  system  (2.1-2),  (2.4-8),  (2.10)  with  a  =  0.001,  i?  =  1cm,  L=20cm, 
dev  =  0.3,  T  =  1200^,  a  flow  rate  of  3505cc,  and  chose  p  =2,  5,  10,  15,  lOtorr.  Our 
model  predicts  that  higher  pressure  improve  the  efficiency  of  the  CVD  system  (see  Figure 
10,  left).  Since  the  mass  flow  rate  is  kept  constant,  higher  pressure  with  constant  tempera¬ 
ture  yield  slower  flow,  higher  impingement  rate,  and  slower  diffusion  rate.  The  diffusion 
rate  and  the  impingement  rate  cancel;  thus,  this  has  the  same  effects  as  varying  the  flow 
rate. 


In  order  to  examine  the  effect  of  temperature  at  constant  dwell  time,  we  selected 
p  =  \Qtorr,  a  =  0.002,  R  =  1cm,  L  =  20cm,  dev  =  0.3,  determined  a  flow  rate  to  main¬ 
tain  a  constant  dwell  time  of  O.OS^cc,  and  solved  the  system  (2.1-2),  (2.4-8),  (2.10)  for 
T  =  800,  900,  1000,  1100,  1200,  1300^!".  Constant  reaction  probability  leads  to  small  vari¬ 
ations  in  efficiency  (see  Figure  11,  left).  This  happens  because  of  the  diffusivity  increase 
with  temperature.  Higher  temperatures  lead  to  a  lower  density  and  precursor  flux  entering 
the  reactor;  thus,  to  a  lower  production  rate  (see  Figure  11,  right). 
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In  order  to  examine  the  effect  of  reactor  length  at  constant  dwell  time,  we  chose 
p  =  lOtorr,  a  =  0.02,  R  =  1cm,  w  =  0.05,  dev  =  0.3,  T  =  1200/sf,  a  flow  rate  determined 
to  maintain  a  constant  dwell  time  of  O.OS^cc ,  and  L  =  10,  20,  30,  40,  50,  60cm .  With 
constant  temperature  and  pressure,  results  of  Figure  12  indicate  that  efficiency  remains 
constant  while  the  production  rate  increases  with  L.  We  show  the  axial  velocity  com¬ 
ponent  Vj.  for  L  =  60cm  using  a  20x20  coarse  mesh  Figure  13  (top-left)  and  using  a  mesh 
with  3  levels  of  h-refinement  in  Figure  13  (top-right).  We  show  a  blow-up  of  the  left- 
lower  comer  of  Figure  13  (top-right)  to  show  the  locally  refined  mesh  that  improves  the 
accuracy  of  the  numerical  solution. 

5.  Discussion  of  Results.  We  have  developed  a  mathematical  model  to  analyze  and  simu¬ 
late  CVD  fiber  coating  processes  for  manufacturing  ceramic  composites.  The  model  is 
used  with  a  general-purpose  adaptive  software  system  for  transient  and  steady  partial 
differential  equations  that  is  capable  of  automatic  mesh  generation,  mesh-order  variation, 
and/or  mesh  refinement.  The  viscous  flow  system  for  the  mixture  is  coupled  with  energy 
equations  for  the  mixture  and  fiber,  a  species  equations  combined  with  a  surface  reaction 
model,  and  a  fiber  coating  thickness  equation.  This  model  has  predicted  the  effects  of 
operating  parameters  such  as  temperature,  flow  rate,  pressure,  etc.  on  efficiency  and  pro¬ 
duction  rates  of  the  CVD  process.  This  model  can  be  easily  altered  to  include  other  effects 
such  as  thermal  diffusion,  multiple  species,  multi-step  surface  reactions,  and  cold-wall  fur¬ 
naces  where  the  fiber  is  heated  by,  e.g.,  electrical  current. 

Concentrating  the  precursor  near  the  fiber  and  away  from  the  hot  wall  is  the  only 
strategy  that  substantially  improves  the  efficiency  of  this  hot  wall  CVD  reactor.  Higher 
reaction  probabilities  yield  higher  efficiencies,  especially  when  the  precursor  enters  near 
the  fiber.  Longer  dwell  times,  controlled  by  either  reactor  length  or  flow  rate,  permit  the 
precursor  to  react  and  deposit  on  solid  surfaces  and,  thus,  improve  the  production  rate  rate. 
While  they  reduce  efficiency,  both  higher  flow  rates  and  precursor  concentrations  increase 
the  CVD  process  production  rate. 
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Future  efforts  will  develop  a  more  realistic  surface  reaction  model  to  include  mutli- 
step  and  multi-species  mixtures.  This  will  be  combined  with  more  experimental  work  to 
understand  the  surface  chemistry.  We  will  also  concentrate  our  effort  on  developing  a 
model  for  a  more  practical  situation  of  coating  tows  of  fibers.  The  gas  infiltration  into  the 
tow  becomes  important  to  obtain  uniform  coatings.  Usually,  when  coating  a  tow  of  fibers, 
the  precursor  reacts  on  accessible  surfaces  of  the  tow  and  closes  the  pores  and,  thus, 
prevents  further  infiltration  of  the  precursor  into  the  interior.  We  plan  to  simulate  and 
analyze  pore  closing.  Software  for  three-dimensional  parallel  adaptive  computations  will 
also  be  developed.  By  Combining  experimental  and  state-of-art  computational  adaptive 
software,  we  hope  to  identify  and  verify  optimal  configurations  and  manufacturing 
processes  much  more  rapidly  than  would  be  possible  by  using  either  paradigm  alone. 
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Table  II.  List  of  Symbols  and  Parameters. 


r,  z 

Radial  and  vertical  coordinates,  respectively 

P 

Mixture  density 

P4 

Density  of  AI2O3 

P/ 

Fiber  density 

P 

Pressure 

Vr 

Radial  component  of  the  velocity 

Vz 

Axial  component  of  the  velocity 

T 

Temperature 

g 

Gravitational  acceleration 

K 

Universal  gas  constant 

Boltzmann  constant 

Avogadro  number 

Viscosity  of  species  i 

Mixture  viscosity 

Cpi 

Specific  heat  of  species  i 

Cp 

Specific  heat  of  the  gas  mixture 

^Pf 

Fiber  specific  heat 

k 

Mixture  thermal  diffusivity 

ki 

Thermal  diffusivity  of  species  i 

'6 

M, 

Fiber  thermal  diffusivity 

Precursor  diffusivity  in  the  mixture 

Molecular  weight  of  species  i 

Mole  fraction  of  species  in  the  mixture 

Yi 

Mass  fraction  of  species  i  in  the  mixture 

^in 

Precursor  mass  fraction  at  inlet 

Pin 

Pressure  at  inlet 

T 

■Lm 

Temperature  at  inlet 

T 

Temperature  at  central  part  of  the  reactor  wall 

Vjru 

Axial  velocity  at  inlet 

M 

Mean  molecular  weight  of  the  mixture 

Rf 

Fiber  radius 

K 

Reactor  radius 

Re 

Reynolds  number 

Fr 

Froud  number 

AH 

Surface  reaction  activation  energy 

Percentage  of  precursor  hat  deposits  on  the  fiber 

Po 

Percentage  of  unreacted  precursor 

F. 

Percentage  of  precursor  that  deposits  on  the  wall 
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Figure  1.  Pressure  (upper-left),  axial  velocity  (upper-right),  temperature 
(lower-left),  and  Precursor  mass  fraction  (lower-right)  resulting  from  the 
solution  of  (2.1-9). 
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Figure  2.  Coating  thickness  h{z)  vs.  z  on  fiber  (left),  thickness  growth  rate 
on  wall  (right)  resulting  from  solving  (2.1-9). 


Figure  3.  Precursor  mass  percentage  deposited  on  wall  (dashed),  fiber 
(solid),  and  unreacted  (dash-dot)  vs.  flow  rate  (left).  Coating  thickness  on 
fiber  vs.  flow  rate  (right)  resulting  from  solving  (2.1-2),  (2.4-8),  (2.10). 
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Figure  4.  Precursor  mass  percentage  deposited  on  wall  (dashed),  fiber 
(solid),  and  unreacted  (dash-dot)  vs.  precursor  distribution  (left).  Coating 
thickness  on  fiber  vs.  precursor  distribution  (right)  resulting  from  solving 
(2.1-2),  (2.4-8),  (2.10). 


Figure  5.  Precursor  mass  percentage  deposited  on  wall  (dashed),  fiber 
(solid),  and  unreacted  (dash-dot)  vs.  reactor  length  (left).  Coating  thickness 
on  fiber  vs.  reactor  length  (right)  resulting  from  solving  (2.1-2),  (2.4-8), 
(2.10). 
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Figure  6.  Precursor  mass  percentage  deposited  on  wall  (dashed), 
fiber(solid),  and  unreacted  (dash-dot)  vs.  reaction  probability  a  (left).  Coat¬ 
ing  thickness  on  fiber  vs.  a  (right)  resulting  from  solving  (2.1-2),  (2.4-8), 
(2.10). 
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Figure  7.  Precursor  mass  percentage  deposited  on  wall  (dashed),  fiber 
(solid),  and  unreacted  (dash-dot)  vs.  precursor  mass  fraction  (left).  Coating 
thickness  on  fiber  vs.  precursor  mass  fraction  (right)  resulting  from  solving 
(2.1-2),  (2.4-8),  (2.10). 
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Figure  8.  Precursor  mass  percentage  deposited  on  wall  (dashed),  fiber 
(solid),  and  unreacted  (dash-dot)  vs.  reactor  radius  (left).  Coating  thickness 
on  fiber  vs.  reactor  radius  (right)  resulting  from  solving  (2.1-2),  (2.4-8), 
(2.10). 


Figure  9.  Precursor  mass  percentage  deposited  on  wall  (dashed),  fiber 
(solid),  and  unreacted  (dash-dot)  vs.  temperature  (left).  Coating  thickness  on 
fiber  vs.  temperature  (right)  resulting  from  solving  (2.1-2),  (2.4-8),  (2.10). 
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Figure  10.  Precursor  mass  percentage  deposited  on  wall  (dashed),  fiber 
(solid),  and  unreacted  (dash-dot)  vs.  pressure  (left).  Coating  thickness  on 
fiber  vs.  pressure  (right)  resulting  from  solving  (2.1-2),  (2.4-8),  (2.10). 
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Figure  11.  Precursor  mass  percentage  deposited  on  wall  (dashed),  fiber 
(solid),  and  unreacted  (dash-dot)  vs.  temperature  at  dwell  time  =  O.OSsec 
(left).  Coating  thickness  on  fiber  vs.  temperature  at  dwell  time  =  O.OSsec 
(right)  resulting  from  solving  (2.1-2),  (2.4-8),  (2.10). 
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Figure  12.  Precursor  mass  percentage  deposited  on  wall  (dashed),  fiber 
(solid),  and  unreacted  (dash-dot)  vs.  reactor  length  at  dwell  time  =  O.OSsec 
(left).  Coating  thickness  on  fiber  vs.  reactor  length  at  dwell  time  =  O.OSsec 
(right)  resulting  from  solving  (2.1-2),  (2.4-8),  (2.10). 


Figure  13.  Axial  velocity  using  a  20  x  20-element  graded  coarse  mesh 
(top-left),  three  levels  of  refinement  (top-right),  and  a  blow-up  of  lower-left 
comer  of  top-left  of  this  figure  (bottom)  resulting  from  solving  (2.1-2), 
(2.4-8),  (2.10). 


-  25  - 


